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Optimizing snake locomotion in the plane. II. Large transverse 

friction 

Silas Alben 

Department of Mathematics, University of Michigan, Ann Arbor, MI, USA 

Abstract 

We determine analytically the form of optimal snake locomotion when the coefficient of 
transverse friction is large, the typical regime for biological and robotic snakes. We find that the 



<: 

\^ • optimal snake motion is a retrograde traveling wave, with a wave amplitude that decays as the 



-1/4 power of the coefficient of transverse friction. This result agrees well with our numerical 
computations. 
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I. INTRODUCTION 



In p| we considered the problem of optimizing snake locomotion in the plane numeri- 
cally. The mechanics of snake locomotion has been described previously by biologists and 



engineers, and modeled by applied mathematicians 



n 



fliii. 



As with other terrestrial locomot- 



ing animals [8], a terrestrial snake pushes against the ground with its body, and obtains a 
reaction force from the ground which propels it forward. However, a wide range of possible 
kinematics can be employed, and determining which are most effective (i.e. efficient) in 
different environments, and why, has been a major theme of locomotion studies 9N20|. 

A first approximation to snake locomotion is to consider motions confined to two dimen- 
sions J2|, yl, |6|. Here the reaction force in the plane of snake motion is due to friction, and 
Coulomb friction provides a simple model. Hu and Shelley |3|, |6| emphasized the impor- 
tance of frictional anisotropy in snake locomotion. In particular, they found that when the 
curvature of the snake backbone is prescribed as a sinusoidal traveling wave, high speed 

transverse friction is large compared to the 



and efficiency is obtained when the coefficient o 
coefficient of forward friction. Jing and Alben 2l| found that the same holds for two- and 
three-link snakes. In biological snakes, the ratio of transverse to forward friction is thought 
to be greater than one [3|,l6|, although how much greater is not clear. In [sl, l6[, a ratio of only 
1.7 was measured for anaesthetized snakes. These works noted that active snakes use their 
scales to increase frictional anisotropy, so the ratio in locomoting snakes is potentially much 
larger, and indeed, a ratio of 10 was found to give better agreement between the Coulomb 
friction model and a biological snake [6|. Wheeled snake robots have been employed very 



successfully for locomotion 22|, |23|, and for most wheels the frictional coefficient ratio (with 



forward 
than 10 



riction defined by the rolling resistance coefficient of the wheel) is much greater 



m. 



In this paper, we develop an analytical solution for the optimal locomotion of a snake in 
the plane, when the ratio of transverse to forward friction is large. We make few assumptions 
on the snake kinematics at the outset, but use numerical computations from [l| to provide 
guiding intuition. We find first that the optimal form of the prescribed backbone curvature is 
a traveling wave. We then find that the root-mean-square (RMS) amplitude of the curvature 
should decay as the —1/4 power of the ratio of transverse to forward friction coefficients. 
Surprisingly, any periodic traveling wave motion can achieve the optimal efficiency, subject 



to the aforementioned restriction on its RMS amplitude, and to the requirement that its 
wavelength, normalized by the snake length, tend to zero. In the limit of large transverse 
friction, the power required to move a snake optimally is simply that needed to tow a straight 
snake forward. 



II. MODEL 



2l| , so we only summarize it here. The 



We use the same frictional snake model as [l|, y, |6, 
snake's position is given by X(s,t) = {x{s,t),y{s,t)), a planar curve which is parametrized 
by arc length s and varies with time t. A schematic diagram is shown in figure [1] 



s = 1 




s = 



FIG. 1: Schematic diagram of snake position, parametrized by arc length s (nondimcnsionaUzed by snake 
length), at an instant in time. The tangent angle and unit vectors tangent and normal to the curve at 
a point are labeled. Vectors representing forward, backward and transverse velocities are shown with the 
corresponding friction coefficients /i/, /ib, and fj,f 



The unit vectors tangent and normal to the curve are s and fi respectively. The tangent 
angle and curvature are denoted 9{s,t) and n{s,t), and satisfy dsX = cos 6', dsH = sin6', and 
K = dgO. We consider the problem of prescribing the curvature of the snake as a function of 
time, K{s,t), in order to obtain efficient locomotion. When K,{s,t) is prescribed, the tangent 
angle and position are obtained by integration: 

s 

(1) 

(2) 
(3) 



e{s,t) = eo{t)+ / K{s',t)ds', 

Jo 

x{s,t) = Xo{t) + / cos6{s',t)ds' 
Jo 

y{s,t) = yo{t)+ / sinOis', t)ds'. 
Jo 



The trailing-edge position (xo,yo) and tangent angle 9o are determined by the force and 
torque balance for the snake, i.e. Newton's second law: 

I pdttxds = / f^ds, (4) 

Jo Jo 

pdttyds = / fyds, (5) 

Jo 

L j-L 

pX^-9tiXrfs= / X^-frfs. (6) 

Jo 

Here p is the snake's mass per unit length and L is the snake length. The snake is locally 

inextensible, and p and L are constant in time, f is the force per unit length on the snake 

due to Coulomb friction with the ground J3[: 

f (s, t) = -pgpt (a3c • n) n - p^ {f^fH{dtX ■ s) + p,{l - H{dtX ■ s))) (dtX ■ s) s. (7) 

Here H is the Heaviside function and the hats denote normalized vectors. When ||i9tX|| = 
we define 9tX to be 0. According to (ITj) the snake experiences friction with different 
coefficients for motions in different directions. The frictional coefficients are /i/, pb, and pt 
for motions in the forward (s), backward (— s), and transverse (i.e. normal) directions (±n), 
respectively. In general the snake velocity at a given point has both tangential and normal 
components, and the frictional force density has components acting in each direction. A 
similar decomposition of force into directional components occurs for viscous fluid forces on 



slender bodies 



25|. 



We assume that the snake curvature K{s,t) is a prescribed function of s and t that is 
periodic in t with period T. Many of the motions commonly observed in real snakes are 



a 



essentially periodic in time |3|. We nondimensionalize equations (jl])-® by dividing lengths 
by the snake length L, time by T, and mass by pL. Dividing both sides by pfg we obtain: 

-— r / dttxds = I fxds, (S) 

H9T' Jo Jo ^ ^ 

-J^.l^-y<'^ = lhds, (9) 

— ^ / X^ ■ du^ds = /" X^ ■ Us. (10) 

/W/t/T^ Jo Jo 

In (I8l)- (fl0l) and from now on, all variables are dimensionless. For most of the snake motions 
observed in nature, L/pfgT"^ ^1 [3], which means that the snake's inertia is negligible. 



By setting this parameter to zero we simplify the problem considerably while maintaining a 
good representation of real snakes. fl8l)- flT0l) become: 

b = (fei,fe2,&3)^ = ; h= [ fjs, (11) 

Jo 

h= [ fyds, (12) 

b3= f X^ ■ ids. (13) 

Jo 

In flTT])-flT3l). the dimensionless force f is 



i{s, t) = -i^ (a3c ■ n) n - UidtX . s) + ^(1 - H{dtX ■ s))\ (^3^ • s) s (14) 



The equations (iTTl) -( lT3l) thus involve only two parameters, which are ratios of the friction 
coefficients. From now on, for simplicity, we refer to fit/fif as fit and fib/fJ^f as fib- Without 
loss of generality, we assume /!& > 1. This amounts to defining the backward direction as 
that with the higher of the tangential frictional coefficients, when they are unequal, fit niay 
assume a.y „o„„egative value. The same model was used in fl fl Q. and was found to 
agree well with the motions of biological snakes in p|. 

Given the curvature k(s, t), we solve the three nonlinear equations ( 1TT|) -([T3l) at each time 
t for the three unknowns XQ{t), yo{t) and 6'o(t). Then we obtain the snake's position as a 
function of time by ([I])-(l3]). The distance traveled by the snake's center of mass over one 
period is 

d = \iij ^(^' 1) - ^(^' 0) d^ + (/ ^(^' 1) - ^(^' 0) ^^) • (15) 

The work done by the snake against friction over one period is 

W= f f f(s,t)-(9tX(s,t)dsrft (16) 

Jo Jo 

We define the cost of locomotion as 

-T (-) 

and our objective here is to find K,{s,t) which minimizes t] as fit ^ oo. 

III. LARGE-^i ANALYSIS 

We now analytically determine the optimal snake motion in the limit of large fif Figure 
|2] shows a few numerical results from pL,], which provide some intuition as we develop the 
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FIG. 2: Numerical traveling-wave optima from [l| for ten different /it ^ 1. a) Snake trajectory over one 
period, from one snake optimization run for /i^ = 1 and /it = 30. b) Curvature versus arc length at the 
instant when a curvature maximum crosses the snake midpoint, for ten fit values: 5, 6, 7, 10, 20, 30, 60, 
100, 200, and 300. c) Snake shapes corresponding to the curvatures in (b). 



analytical solution. In [1] we show that for yU^ ^ 6, the numerically-computed optimal 
motions are always (retrograde) traveling waves, an example of which is shown in panel a 
for fit = 30. A sequence of snapshots of the snake is shown over a period of motion. The 
snake moves from left to right, and appears to follow a sinusoidal path, very similar to 
what has been observed biologically and studied numerically |6] . We have computed similar 
traveling- wave optima at a range of /i^ ^ 1. In panel b we compare the curvatures from 
these optimal motions at a particular instant when a curvature maximum occurs at the 
snake midpoint, for ten different values of fif 5, 6, 7, 10, 20, 30, 60, 100, 200, and 300. The 
profiles are similar, but the amplitudes decay monotonically as fit increases. Panel c shows 
the ten snake shapes corresponding to the curvatures in panel b. Consistently, we see the 
decay of the traveling- wave amplitude as fit increases. With this picture from the numerics, 
we now proceed to derive the analyical solution in the asymptotic limit of large fit- 
We begin with the position of the snake in terms of its components: 



X{s,t) = {x{s,t),y{s,t)). 



(18) 



If fit is large, the snake moves more easily in the tangential direction than in the transverse 



direction. Furthermore, the most efficient curvature function will give motion mainly in the 
tangential direction, to avoid large work done against friction. We assume that the mean 
direction of motion is aligned with the x-axis, and that the deflections from the x-axis are 
small — that is, \y\, \dty\., \dsy\ and higher derivatives are 0{fi^) for some negative a. This has 
been observed in the numerical solutions, and makes intuitive sense. Small deflections allow 
the snake to move along a nearly straight path, which is more efficient than a more curved 
path, which requires more tangential motion (and work done against tangential friction) for 
a given forward motion, i.e. d. 

We expand each of the terms in the force and torque balance equations in powers of \y\ 
and retain only the terms which are dominant at large fif We expand the cost of locomotion 
similarly, and find the dynamics which minimize it, in terms oi y{s,t). 

We ffist expand x{s,t). We decompose x{s,t) into its s-average x and a zero- s- average 
remainder: 

.(M)^^(M)+/'cos9(/,«)..' (19) 

where -f means the constant of integration is chosen so that the integrated function has 
zero s- average. Hence 

dtx{s,t) = dtx{s,t) + J- -dte{s',t)sme{s',t)ds'. (20) 

We denote the s-averaged horizontal velocity (the horizontal velocity of the snake center of 
mass) by 



U{t) = dtx{s,t). (21) 

We expand the integrand in (1201) in powers oi \y\ and its derivatives: 

- dtOis, t) sin e{s, t) = -dsty{s, t)dsy{s, t) + 0{\dsy\^). (22) 

The quartic remainder term in (I22p . 0{\dsy\^)^ actually includes other quartic terms involv- 
ing time derivatives of dsy, but for brevity we use the assumption (which will be correct for 
our solution) that \y\ and all its derivatives are of the same order of smallness with respect 
to Hf We define 

h{s,t) = -[ -ds'ty{s\t)dMs\t)ds' (23) 

and with the small-defiection assumption, (!20|) becomes 

dtx{s,t) = U{t) + h{s,t) + 0{\dM''), (24) 



and 



dtX{s, t) = {U{t) + h{s, t), dty{s, t)) + 0{\dsy\ 



(25) 



We will see that for a general class of shape dynamics near the optimum, U{t) is 0(1), i.e. 
of the same order as one snake length per actuation period, even in the limit of large /ij. U 
is therefore the dominant term in f l24l) and ( 125|) . We use this assumption on U to expand 
the velocity 2-norm: 



\dtX{s, t) II = ^dtx^ + dty^ = uU + ^ + ]^^\ + O 



M')- 



(26) 



In writing 0(|y|^) we have again lumped terms with four or more powers oi \y\ and/or its 
s- and t-derivatives. Dividing ( 125|) by ( !26|) . we obtain the expansion for the normalized 
velocity: 



/ 



atX(s,t) 



i-^T?+o(i,n 



\ 



\ U \ f/ 2 t/2 
We next expand the unit tangent and normal vectors: 



0{\y\ 



(27) 





n 



1 - kdsy' 
dsy 

-dsy 
1 - \dsy' 




(28) 
(29) 



Using ([27)) - ([29 



1 



dtV 



5*X.n= --9., 1---^ 



^ty f ^ '^a 2\ , nn \ 



(30) 

(31) 



Using f l30|) and fl3Tl) we can write the force and torque balance equations flTT]) -f fT3|) and 
the cost of locomotion rj in terms of U, h, and derivatives of y. 

A. Traveling-wave optimum 

We first expand rj and use it to argue that the optimal shape dynamics is approximately 
a traveling wave. This allows us to neglect certain terms, which simplifies the formulae for 



8 




the rest of the analysis. If the power expended at time t is P{t), the cost of locomotion is 
ri= [ P{t)dt I [ U{t)dt. (32) 

^^rfsdt / / U{t)dt. (33) 

^(^) / Jo 

In ( l33l) we have assumed that the tangential motion of the snake is entirely forward, with 
no portion moving backward, so the coefficient in front of the tangential term is unity (i.e. 
fif, rather than a term involving both fif and fib). This assumption is not required but 
it simplifies the equations from this point onward and will be seen to be correct shortly. 
The numerator of (133|) can be separated into a part involving tangential velocity and a part 
involving normal velocity with a factor of fit- Using ( l26l) . ( 130|) . and ( 13T]) . we can list the 
terms with the lowest powers of y in each part: 

^■s)'^ = l + 0(|yp) (34) 

fit {dtX-hy i^ = fit (d,y - ^) + Oifitlvn (35) 

The optimal shape dynamics is a y{s, t) that minimizes t] subject to the constraints of force 
and torque balance. The leading-order term in flM]) is 1, which is independent of the shape 
dynamics. At the next (quadratic) order in y, terms appear in both (^^ and fl35|) : only that 
in (!35|) is given explicitly. It is multiplied by fit, so it appears to be dominant in the large 
fit limit. Thus, to a first approximation, minimizing rj is achieved by setting 

dsy - tt) = 0, (36) 

so a first guess for the solution is a traveling wave, i.e. a function oi s + Ut with U constant. 
However, such a function does not satisfy the x-component of the force balance equation 
(ITTj). which is: 



n(^dtX-s^s, + fit{dtX-njn.,^ds = 0. (37) 

Inserting (!30|) and (!3T]) and retaining only the lowest powers in y from each expression, 

J 1 + fit (dsy - ^) dsv ds = 0. (38) 

which cannot be satistied by a function of s + Ut. Physically, the first term (unity) in ( 138|) 
represents the leading-order drag force on the body due to forward friction. For a function of 



s + Ut, in which the deflection wave speed equals the forward speed, the snake body moves 
purely tangentially (up to this order of expansion in y) along a fixed path on the ground, 
with no transverse forces to balance the drag from tangential frictional forces. Also, when a 
function of s + Ut is inserted into fl5Sl) all factors of U cancel out and U is left undetermined. 
This is related to a more fundamental reason why posing the shape dynamics as a function 
oi s + Ut is invalid: the problem of solving for the snake dynamics given in section [TT] is not 
then well-posed. In the well-posed problem we provide the snake curvature and curvature 
velocity versus time as inputs, and solve the force and torque balance equations to obtain 
the average translational and rotational velocities as outputs. Thus U{t) is one of the three 
outputs (or unknowns) that is needed to solve the three equations. 
Instead we pose the shape dynamics as 

y{s,t)=g{s + Ua), (39) 

which is a traveling wave with a prescribed wave speed 1/^, different from U in general. 
Since y is periodic in t with period 1, g is periodic with a period of Uw Inserting fl39|) into 
(138|) . we obtain an equation for U in terms of f/^ and g: 



i + ^it(i-jf] [ g\s + U^tfds = 0. 



(40) 



In fjlU]) . as yUi — )■ oo, (1 — U^/U) — > 0". In the solution given by fl5I?]) and (HUl) . the shape 
wave moves backwards along the snake at speed Uu,, which propels the snake forward at a 
speed [/, with U slightly less than Uyj. Therefore the snake slips transversely to itself, in 
the backward direction, which provides a forward thrust to balance the backward drag due 
to forward tangential friction. We refer to (1 — Uw/U) in (140!) . a measure of the amount 
of backward slipping of the snake, as the "slip." In figure [2|i, the snake slips transversely 
to itself, and backward (leftward), even as it moves rightward. Here the slip is fairly small 
because /if is 30, fairly large. 

We now solve fj40|) for (1 — U^/U) and insert the result into fl33|) to find g and U^ which 
minimize rj. Using the lowest order expansions (IMjl and (!35l) with (HOl) we obtain: 



10 



where 



rj = I P{t)dt 



U{t)dt. 



\{t) f\ + ^,{d,y-^] +0{\y\\^xt\y\')dsdt 



U 



U{t) dt. 



-dt. 



(41) 
(42) 

(43) 



1 + 



l^t{9'{s + Ua)^) + 0{\g\\i^t\9n 



(44) 



{g'is + U^tf) = / g'{s + U^t)^ds. 
Jo 

If we minimize (H3!) over possible g for a given /it, using only the terms in the expansion that 
are given explicitly, we find that r] tends to a minimum of 1 as the amplitude of g diverges. 
However, in this limit our small-y power series expansions are no longer valid. We therefore 
add the next-order terms in our expansions of t] and the Fx equation and search again for 
an r^-minimizer. 



B. Optimal wave amplitude 



We now expand to higher order the components of rj, ( IMl) and (l35l) 

(45) 
(l + Odt/l^)), (46) 



dtX-s] ^^ = l-Uy-jf] + Jj + 217^ + Oi\y\') 



Ht dtX ■ n 



Uit) 



Uit) 



/ii 



9ty a ^\ , dty f h 1^ 2 



and equation flTTj) . 



^-5.y +-^ -^-2^.y 



dty , <9iy Z' /i , 1^2 



1 + /i* ( 5.1/ - -^ + -^ ( ^ + 2^^?/ ) ) dsV ds = 0. 



(47) 



We again assume the traveling wave form of y (!39|) , and use this to evaluate h given by ( !23|l 
in terms of (7: 

/i(s, t) = -^g'{s + f/^t)2 + ^{g'{s + t/^t)^). (48) 

Inserting the traveling-wave forms of y ( !39l) and /i ( l48l) into ( l47l) we obtain 



l + /it 



l-^)(^'(^ + t/.i)') + ^(^'(^ + f/.t)2)' 



0. 



(49) 



11 



To obtain ( H9|) we've used the fact that the shp (1 — Uw/U) — )■ as /it — )■ oo, which will be 
verified subsequently. One can also proceed without this assumption, at the expense of a 
lengthier version of (1491) . We insert the traveling- wave forms of y and h into (1451) and (H^ 
to obtain rj, again with the assumption that (1 — Uw/U) — )■ as /x^ — )■ oo to simplify the 
result: 

V = ^ + l' lig'is + Ua?) +fit(^l-^ + ^{g'{s + U^ty)\ {g'{s + U^tf) dt (50) 

Solving for the slip (1 — Uw/U) from (I49p and inserting into (I50p . we obtain an improved 
version of (HHl): 




r] = l / dt. (51) 



^ + ^<«''- + ^-'» %.M. + c/„y) ^""^'^'^''^''' 



( 15T1) is a more accurate version of ( l43l) . with an additional term in the denominator of 
the integrand which penalizes large amplitudes for g. If we approximate {g'{s + U^tY) as 
constant in time, we obtain 

^ = 1 + \{9") + ^w^ + 0{\g\\^it\g\% (52) 

which is minimized for 

{g'Y' = 2VV,-^/^ (53) 

Thus the RMS amplitude of the optimal g should scale as /i^ . fl53|) represents a balance 
between competing effects. At smaller amplitudes, the forward component of normal friction 
is not large enough to balance the forward component of tangential friction, so the snake 
slips more, which reduces its forward motion and does more work in the normal direction. 
At larger amplitudes, as already noted, the snake's path is more curved, which requires more 
work done against tangential friction for a given forward distance traveled. (|53|) is only a 
criterion for the amplitude of the motion. To obtain information about the shape of the 
optimal g, we now consider the balances of the y-component of the force and the torque. 

C. Small- wavelength shapes 

So far we have searched for the optimal shape dynamics in terms of y{s, t), but in fact it is 
the curvature k{s, t) which is prescribed. We obtain y{s, t) from the curvature by integrating 

12 



twice in s: 



y{s,t) = y{0,t) + sm9{s',t)ds' 
Jo 

= y{0,t)+ f e{s',t)ds' + 0{y') 



y{0,t) + 



6(0, t)+ / K{s",t)ds 




s rs' 



ds' + 0{y^). 



= y{0,t) + e{0,t)s+ / K{s",t)ds"ds' + 0{y''). 
Jo Jo 

= Y{t) + sR{t) + k{s,t) + Oiy""). 



(54) 
(55) 

(56) 

(57) 
(58) 



where Y and R are defined for notational convenience. Prescribing tlie curvature is equiva- 
lent to prescribing k{s,t). We set 



k{s,t) = g{s + Uu,t) 



(59) 



so we liave tfie same form for y as before, with an additional translation Y(t) and rotation 
R(t). Y and R are determined by the y-iorce and torque balance equations: 

dtX -s) Sy + fiti dtX ■n]ny 



dt'X ■ sj {xsy - ys^) + ^t ( f^tX ■ nj {xuy - yn.. 
We expand these to leading order in y and obtain 



ds = Q. 
ds = 0. 



(60) 
(61) 







U U \U 2 
dty , dty f h ^ 1^2 



^■^[9sy-i^ + i^[jj^-dsyn]ds^o. 



(62) 
(63) 



We insert ( 158|) with (!59|) into the three equations (1471) . (l62l) . and (163|) to solve for the three 
unknowns, U, Y, and R in terms of g and t/^. We obtain: 



1+f^t 



f/. 



1 



-- + l + ^(.'^)).'^ 






rfs = 0. 



' '-% + i + ^(^'')l^'-^-^ + i?rf3 = o. 



u 



u u 






ds = 0. 



(64) 
(65) 
(66) 
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We solve (165|) and (166!) for Y and R in terms of U: 



Y' 

— -R = AB-6A (67) 

/?' 

— = 12A - 65 (68) 






(69) 
(70) 



We then solve (164|) for [/ in terms of (7: 



r/ 1 1 



r^ = 1 + ± ((^'2^ ^ /y^X 

We now form 77, ( 133|) with terms given by (145|) and ( 146|) . using the updated from of y (( 158|) 
with f l59|) ) and f lTT]) . We obtain a more correct version of f l5T|) : 

ri=l/ (it. (72) 







We can find ?7-minimizing gf in a few steps. First, let {L^} be the family of orthonormal 
polynomials with unit weight on [0, 1] (essentially the Legendre polynomials): 



(73) 



f LiLj ds = 6ij ; Lo = l, Li = Vl2{s - 1/2) , . . . . 
Jo 

At a fixed time t we expand g' in the basis of the L^: 

CXD 

g\s + U^t) = Y,Ck{t)Lk{s). (74) 



fc=0 



Then we have: 



f) = J2'k{tr, (75) 



fc=0 
7\2 _ ^ /'j-\2 



(^r = Co(t)^ (76) 

3i{g')~2{sg')f = {L,gy = c,itf. (77) 

Inserting into fl72]) . rj becomes 

7] = !/ dt. 

(78) 
14 



The denominator of the integrand in (178!) is minimized when 

oo 

Co(t) = ; ci(t) = ; ^Cfc(t)2 = V2^;^'\ (79) 

k=2 

The function g{s + U^t) that minimizes rj is that for which fl79|) holds for all t. Then the 
integral in ( !78|l is maximized, so r^ is minimized. Recall that g is a periodic function with 
period Uy^. The relations (17^ hold for any such g in the limit that f/^ — t- 0, as long as g is 
normalized appropriately. Define 



1 



U^ \ V2 



A^[^J^ g'{x)'dxj . (80) 

Then 

/ g'{s + Ujfds = A^ + 0{U^) (81) 

Jo 

[ g\s + U^t)ds = 0{U^) (82) 

/ sg'is + U^t)ds = 0{U^). (83) 

Jo 

( IHTj) - ( l83l) are fairly straightforward to show, and we show them in Appendix [A] for com- 
pleteness. If (ED) - (IMl) hold, then in the limit that f/^ ^ 0, ([79]) holds with A = 2^/Vt~^^^ 
by (!75|) - (177|) . A physical interpretation of the small- wavelength limit is that the net vertical 
force and torque on the snake due to the traveling wave alone (g) become zero in this limit, 
so no additional heaving motion (Y) or rotation (R) are needed, and thus the additional 
work associated with these motions is avoided. Similar "recoil" conditions were proposed 
as kinematic constraints in the context of fish swimming |15i]. We also note two other scal- 
ing laws. For a given gf ~ ^u^ , by (I7T]) the slip (1 — Uyj/U) ^ ^t ^'^^ by (!67j) - (!70ll . 

D. Optimal cost of locomotion and comparison with numerics 

To obtain an 77- minimizing function g*, let g be any periodic function with period Uw 
Define 

g\s + U^t) = ^ g{s + U^t) (84) 

where A is given in ( IHOj) . Then in the limit f/y, — )■ 0, 

r/(/)^l + v^/xr'/VO(/x,-i), (85) 
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by the estimate in ( !72|) . Thus the optimal travehng wave motions become more efficient 
as jjit — )■ oo, which agrees with the numerical solutions in [ij. The numerically-determined 
optima in figure [2] do not have small wavelengths, due to the finite number of modes used. 
This is explained further in the Supplementary Material of [l|. 
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FIG. 3: The /it-scalings of the numerical traveUng-wave optima from figure [5 and [l|. (a) The cost of 
locomotion 77 minus unity, on a log scale, for the numerical optima (squares) of [iJ together with rj — 1 for 
the analytical optimum (solid line) . (b) The curvature from figure ^}p rescaled by ^^ . (c) The deflection 
from figure [H; rescaled by /x/ . 



We compare the numerically-found optima to the theoretical results in figure [31 Here 
we take the plots of figure [2] and transform them according to the theory. In figure ^ we 
plot r] — 1 versus fxt for the numerical solutions (squares) from 'l|], along with the small-t/^ 
theoretical result (solid line), which is given by ( l85l) . The numerics follow the same scaling 
as the theory, but with a consistent upward shift. Some degree of upward shift is expected 
from the finite-mode truncation in the numerical computation, which makes the numerical 
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result underperforin the analytical result, as explained in [l|. The circle shows the efficiency 
for a curvature function with a higher wavelength than those which could be represented in 
the numerical optimization, and its distance from the solid line is 0{fi^^), of the order of 
the next term in the asymptotic expansion of r/ — 1. 

Figure [Sb shows the curvature, rescaled by /i/ to give a collapse according to (1551) . The 
data collapse well compared to the unsealed data in figure [2)d. Figure [St shows the body 
shapes from figure [2]: with the vertical coordinate rescaled, and the shapes plotted with all 
centers of mass located at the origin. We again find a good collapse, consistent with the 
curvature collapse. 

IV. CONCLUSION 

We have studied the optimization of planar snake motions for efficiency, using a fairly 
simple model for the motion of snakes using friction. When the coefficient of transverse 
friction is large, our analysis shows that a traveling-wave motion of small amplitude is 
optimal. The amplitude tends to zero as the transverse friction coefficient tends to infinity, 
scaling as the transverse friction coefficient to the -1/4 power, and the efficiency tends to 
unity in this limit. The corresponding power is that for a straight snake towed forward. 

In l|] we were able to compute many optimal motions at moderate and small fit also. 
Those found at small fit also corresponded to traveling waves (direct now, not retrograde), 
and at zero /it, a simple triangular- wave solution was found with optimal efficiency. It 
remains to determine optimal motions with small but nonzero fit analytically. At moderate 
fit, a region of standing-wave or ratcheting motions was found. It may be possible to compute 
some of these motions analytically using a set of assumptions specialized to the moderate- /xj 
regime. 
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Appendix A: Small-^vavelength integrals 

To show (ISTl) - (183|) . we first note that since g is periodic, g' is periodic with zero average 
over a period. To show flHTl) . we decompose the integral in fIST]) into two intervals: 

/ ^'(s + U^tfds = / c/'(s + U^tfds + / ^'(s + U^tfds. (Al) 

The first integral on the right of (lAip is over an integral number of periods of g'"^, and is 
thus [l/Uu]\Uu,A^ = A? + 0{Uu))- The second integral is over an interval of length < f/^, 
and is thus O ([/„,), which shows (15T1) . To show (152]) . we use (I AID again but with ^f' in place 
oig": 

/ g'{s + Uu,t)ds= / ^'(s + ?7^t)rfs + / g'{s + U^t)ds. (A2) 

The first integral on the right of ( ]A2p is over an integral number of periods of g', and is thus 
zero. The second integral is over an interval of length < [/^, and is thus 0{Uw). To show 
(183|) . we decompose the integral into three parts: 

/ sg'{s + Uwt)ds = / f/u;[s/t/u,j5f'(s + f/i„t)(is 

Jo Jo 

+ / {s - U^[s/U^\)g'{s + U^t)ds + / V(s + t/«,t)rfs. (A3) 

The first integral on the right of flA3p substitutes a constant approximation to s on each 
subinterval of period length. It is identically zero since g' has mean zero. The second integral 
is the error in the constant approximation to s, 0{Uyj) on each subinterval, of which there 
are 0(l/f/^), each of length [/^. Thus the second integral is 0{Uu,)- The third integral is 
also 0(t/u,), as explained for (jAip and (]A2p . so (!83|) holds. 
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